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Abstract 



Rho photo-production on complex nuclei is re-examined using a gener- 
alised vector dominance model which succesfully predicts the oberved nuclear 
shadowing in real photoabsorption and deep inelastic scattering. This model 
is shown to give a good fit to rho photoproduction data on both nucleons 
and complex nuclei, in which the disagreement between the measured 7 — p 
coupling and the 7 — p coupling required by the simple vector dominance 
model is eliminated. The pN total cross-sections required are similar to those 
predicted by the additive quark model; and the magnitude of the correction 
to simple vector dominance is consistent with that inferred from the analysis 
of real photoabsorption and deep inelastic scattering. 
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I. INTRODUCTION 



In this paper, we address the apparent contradiction between the well-known simple vec- 
tor dominance (SVD) treatment of rho photoproduction; and the very succesful generalized 
vector dominance (GVD) treatment of nuclear shadowing in real photoabsorption and deep 
inelastic scattering. 

In 1989, one of us M pointed out that the observed qualitative features of nuclear shad- 
owing in deep inelastic scattering [H are simply and naturally accounted for in a generalized 
vector dominance(GVD) model || which can be shown to be dual to the parton model [|T],|J. 
The crucial feature of this model, in the context of shadowing, is that the cross-sections cvn 
for scattering a sequence of hadronic vector states V = p,p',... from nucleons are required 
to be approximately independent of their mass my. This then leads to approximate scal- 
ing behaviour for shadowing and a rapid decrease in the effect as x increases from zero, as 
observed in the data 0. Somewhat later, in 1993, the same model was shown ||[] to give 
an accurate quantitative account of both real photoabsorption on nuclei and the precise 
shadowing data that had by then become available for virtual photons || . 

The above GVD model is consistent with the fundamental QCD picture of strong in- 
teractions in appropriate kinematic regions^. It has implications for rho photoproduction, 
because it necessarily includes substantial contributions from "non-diagonal" diffraction dis- 
sociation processes of the type VN — > V'N in addition to "diagonal" elastic processes of 



1 This paper was completed while the standard review |Q| of nuclear effects in structure functions 
was in press. It is therefore not covered in this review, despite the fact that it was published slightly 
earlier. 

2 This is discussed explicitly in recent papers || in which GVD models are extended to incorpo- 
rate the small- x rise in the proton structure function associated with the "hard Pomeron", which 
is observed at large Q 2 at HERA. The resulting predictions for its behaviour at low Q 2 have 
subsequently been confirmed by experiment [|l0|]. 
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the type VN — > VN. This feature is also implied by duality with the parton model In 
addition, without such processes, the cross-sections ayN would be forced to decrease rapidly 
with mass my to maintain approximate scaling in the nucleon structure functions; and shad- 
owing would die away as Q 2 increases at fixed x, in contradiction to approximate scaling 
for the nuclear structure functions. These diffraction dissociation processes also contribute 
to r/io-photoproduction, as illustrated in Figure [l], giving significant corrections to the well- 
known predictions of the simple vector dominance(SVD) model. This result is potentially 
a serious problem, since rho photoproduction on both nucleons and nuclei has long been 
regarded as the outstanding success of the SVD model; and if the agreement between the 
SVD predictions and the data were really good, it would clearly undermine the above GVD 
approach. However, as Donnachie and Landshoff ]lTJ have recently pointed out, the SVD 



predictions for the cross-sections on nucleons are lie approximately 16% above the measured 
values. 

In this paper we shall investigate rho photoproduction on both nucleons and nuclei within 
the framework of the GVD model used to successfully account for shadowing in real and 
virtual photoabsorption. The aims are to see whether model can resolve the discrepancy 
between the predictions of SVD and the nucleon data, while retaining the successful predic- 
tions of SVD for the A dependence on nuclei; and if so, whether sign and magnitude of the 
diffraction dissociation terms required by the rho photoproduction data are consistent with 
those required by the real photoabsorption and structure function data. 

II. RHO PHOTOPRODUCTION ON NUCLEONS 

We begin by reviewing rho photoproduction on nucleons 

7 + N -> p + N (N = n, p) 



-a topic in which there is renewed interest because of recent data from HERA |l2|,|i"3|j . At 



high energies, the process is usually described in terms of the SVD model of Figure [I], leading 
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to the well-known relationP] 

/ 7 ,(M = 0) = f/ w (s,t = 0) (SVD) (1) 

Jp 

where f a b is the scattering amplitude for a + N — > b + N and f p is the 7 — p coupling. Using 
the optical theorem, this gives 

da , „ s 47r 



^ 7 ( S) t = 0) = «^a p 2 iV [l + r / 2 ] (SVD) (2) 



for the forward differential cross-section, where <j p n is the total cross-section for p N scatter- 
ing and 77 is the ratio of the real to imaginary part of the forward p N scattering amplitude. 
The value of <j p n is usually taken from the additive quark model prediction 

VpN = ^ \ a TT+N + &it- n] ; (3) 

and 77 is estimated from Regge pole ideas. 

The resulting SVD prediction (0, ^) was long ago compared with experimental photo- 
production data on protons available at energies below 20 GeV to yield a value |TJ] 



^ = 2.44 ±0.12. (4) 

47T 

for the coupling constant. At the time, this simple picture was consistent with the cor- 
responding SVD analysis of vector meson photoproduction on nuclei |15|JT6[] and with the 
then not very precisely known value of the 7 — p coupling obtained directly from the mea- 
sured decay width Y(p — > e + e~). Since then events have moved on and it is important to 
check whether this consistency still obtains. The precision of the measured decay width 



T(p — > e + e ) has greatly increased |L7j, and now gives 



f 2 

J o 



j*- = 2.01 ±0.10 , (5) 

4 7T 



3 For an review of vector meson photoproduction and other diffractive photoprocesses in the con- 
text of SVD, see Leith 



which is not in good agreement with the phenomenological value (|]). This problem is 
confirmed by a recent comparison |IT| incorporating higher energy photoproduction data 
fl2| , [T8"|| , which finds that the SVD prediction (0, |3[) obtained using the directly measured 
coupling (^) lies on average about 16 % above the measured data. 

The discrepancy between the experimental data and the SVD approximation presum- 
ably arises from contributions which are neglected in this approximation. There are two 
obvious possibilities. The first is finite width effects associated with the p° — > 7r + 7r~ decay 
channel, as shown in Figure |2|. Detailed calculations show that this is very unlikely, since 
the correction arising from the modification to the rho propagator is largely cancelled by 
corrections arising from the pion scattering terms |L9| . 

The second possibility is to attribute the discrepancy to neglected contributions from 
higher mass vector states, leading to the generalized vector dominance (GVD) model shown 
in Figure |l[ In principle many states could contribute, but in practice only the lightest states 
are expected to be important, since the diffraction dissociation amplitudes f(p'N — > pN), 
and to a lesser degree the couplings (e/ jy), are expected to decrease rapidly with increasing 
2D|. Hence it is reasonable to approximate these contributions with that of a 



mass rripi 



single effective p' state, when the SVD prediction ([]]) is replaced by 

J P Jp' 

This form also follows directly from the generalised vector dominance model || used to 
successfully predict the observed shadowing effects in both real photoabsorption and deep 
inelastic scattering on nuclei ||. We shall not discuss this model further, but just use it to 
specify the properties of the effective p', which are 

m 2 pl = 3m 2 p fff = m p ,f p /m p (7) 

together with 

fp'p' = fpp fp'p = fpp' = —tfpp (8) 

for the forward scattering amplitudes, where 



e « 0.2 . (9) 

Here, we shall retain (|6] - ^) but treat e as a free parameter to be determine by the p- 

photoproduction data. The forward cross-cross-section is then given by 

2 

(GVD). (10) 

For fixed a P N and r/, the results of SVD with the effective coupling (§]) are reproduced by 
GVD with the physical coupling @ and e = 0.160 ±0.055, which is in good agreement with 
the value @ quoted above. 

We thus have two accounts of the data on p photoproduction on protons: either the 
SVD model in which the 7 — p coupling is adjused to fit the data; or the GVD model in 
which this coupling is fixed at its measured value and and the parameter e is adjusted to 
fit the data. The latter is more compatible with the ideas used to understand shadowing in 
nucleon structure functions, but its predictions for p photoproduction on nuclei have never 
been examined. It is our aim to remedy this omission and discuss the implications of the 
results. 



d o . 
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III. RHO PHOTOPRODUCTION ON NUCLEI 



In the generalised vector dominance model, an incident photon may convert into a a whole 
sequence of isovector, vector mesons V — p, p', p", . . . while traversing a nucleus. According 
to Glauber theory [^l|] [p2], this possibility is dealt with by introducing appropriate optical 
potentials U^y, ^Vp Uvv, Uyv'i which characterise the features of the scatterer, where V, V 
are arbitrary members of the vector meson sequence. The resulting wave equation then 
reads 



\ frr ut\ 



V 







V 



£/ 7 U 







(ii) 



where 
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{V)yy, =5 VV/ (V 2 + k 2 V ) (12) 
\U)yy, = UyV 

(*)v =Vv 

In this paper we restrict ourselves to two hadronic channels p, p' and use the eikonal 
approximation 

^ 7 = $ 7 e ifc ^ # p , p , = Q^e*".'* (13) 

where the reduced wave functions $ 7 , $ p , $ p / are assumed to vary slowly enough for all 
second derivatives to be discarded. Furthermore, since we are only interested in the photo- 
production amplitude, it is sufficient to work to order 0{\fa) only. Hence, when the photon 
wavefunction occurs multiplied by factors U yp or U 77 , which are already of this or higher 
order, we can replace it by the incident photon wave \I/ 7 = e tkjZ . With these approximations, 
the reduced wavefunctions of the vector mesons satisfy the coupled differential equations 

jU p (b, z) = ~ (u ip e^ + U pp $ p (b, z) + EV*„(b, z)e V) (14) 

^$ p ,(b, z) = (U^jV + C/ pV ^(b, z) + U p/p %(b, z)e n ^ z ) (15) 

where the space coordinates are parameterized through the z-coordinate in the direction of 
the incident photon wave and the impact parameter b in the plane perpendicular to the 
z-axis. The expressions q\\ i:j = k\ — kj denote the minimal longitudinal momentum transfer 
between particle i and j, where i and j are members of the set 7, p, p'. 

A. Optical potentials and nuclear densities. 

The optical potentials are given by [^]J 

Uij = -Aixfij n(b, z) Uy = Uji 
7 



(16) 



where the forward scattering amplitudes on nucleons are given by eqs.(|] - ||) together 
with the GVD analogue of (|6|) for the p': 

f-yp 1 = T$pp' + Tfp'p' ■ ( 17 ) 
Jp Jp' 



Two different models JEJ are employed for the nuclear density n(b,z), depending on the 
size of the nucleus. 

• A > 16 

For heavy nuclei, A > 16, we use a Fermi gas-like mass distribution of the form 

n{r) =n w (l + exp ( LZ ^ L )) ' ■ (18) 

where the 'skin-thickness' c = 0.545 fm. The Woods-Saxon radius Rw is given as the 
solution of the equation 

A f(R w ) = n w = A Ph f(R Ph ) (19) 
corresponding to a fixed central density riw, where 

M?) ( 1 + 7r 7 

and Rpb = 6.626 fm. 



2\ - 1 
,cr N 



M=\7Z3i ! + ^ (2°) 



A < 16 

For light nuclei, the Woods-Saxon formula is not a good description. Instead we use a 
shell model (harmonic oscillator) density, given by 



where 



„ (r) = n , |1 + , ( _ j exp _(_) , (21) 



S = ^ n. = (0?fl S )- 3 (22) 



and the shell model radius i?,<? = 0.708 A» fm. 



Finally, to take two-body correlations into account, we modify the nuclear density functions 
by the replacement |23] 



mr) _^ n ( r ) | i + I i c apN n ( r ) ^L^T j | (23) 
where the two body correlation length l c =0.3 fm. 



B. Evaluation of the cross-sections 



To evaluate the cross-section for p photoproduction on heavy nuclei, we need to solve 



([14], |15|) with the initial conditions 

® p (b,z = -oo) = 5y(b,z = -oo) = (24) 
corresponding to incident photons. The forward differential cross section is then given by 

^(* = 0) = ^I^p(0)| 2 . (25) 
where the forward scattering amplitude F 1P is related to the "profile function" 

r p (b) = -lim<l>,0M) . (26) 



by the standard result [21 



poo 

F lp (0)=ik p dbbT(b). (27) 

J 

From ([14], [15]), it is clear that the reduced wavefunctions only vary with z at fixed 
impact parameter b in regions where optical potentials (|i~6|) are non-zero. Since the nuclear 
densities fall off very rapidly beyond the nuclear radii, we neglect such variations for \z\ > 
4:Ryy, measured from the centre of the nucleus. The initial conditions ( [24]) are imposed 
at z — —4Rw and (0, |15|) are integrated up to z = +4Rw at fixed impact parameter b 
using a fourth-order Runge-Kutta algorithm with an adaptive stepsize control. The profile 
function is then evaluated at z — +ARy/ rather than infinity and used to compute cross- 
sections, which are compared with experiment in the concluding section. Typical results for 
the reduced wavefunction <3> p and the profile function T(b) are shown in Fig. ^| and Fig. |] 
respectively, verifying the assumed disappearance of structure by 4Rw- 



C. A simple approximation 



Before comparing with experiment, we comment briefly on an approximation that gives 



a useful check on our numerical solutions. In the SVD model, the eikonal equation (14) 
with f pp > = can be integrated explicitly and this result is easily generalized to the coupled 
channel GVD case (|TJ|, |T^) in the high energy limit when 

e""'-' z = 1 . (28) 

This approximation has been used to analyse the p photoproduction data by Kroker 



giving, for example, a value <j p n = 29.0 ±1.2 mb corresponding to an incident photon energy 
of 6.1 GeV. The corresponding prediction for the nuclear cross-sections is shown in Figure |5], 
where it is compared both to the data and to our "exact" predictions for the same parameter 
values. In fact 

e iq K P ' z = i for £ 7 = 6.1 GeV, z « 3 / . 
Since z = 3 f is well within a large nucleus for small impact parameters b ~ 0, it is clear that 



(Pq) is quantitatively unreliable at the energies where data currently exists. We shall not 



consider it further, except to note that our numerical solutions reproduce both the standard 



SVD predictions and Kroker's approximate results |24j when we impose f pp i = or fl28|) 
respectively. 

IV. RESULTS AND CONCLUSIONS 

In this section we present the results of an optical model analysis of the experimental 
data on complex nuclei using the GVD model described. At each energy the cross-sections 
depend on four parameters: 

P 

Our strategy is to fix the coupling f p at its measured value @ and the phase ratio rj at the 
values implied by the Regge pole parameterization 
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— I — p-naip —1 — p~' K0L R 

f pp (s, t = 0)= g P — -s ap + g R — -s a « (29) 

sm{7cap) sm{7caji) 

where the subscripts P and R denote the pomeron and Regge contributions respectively and 
the constants op, Or, gp and gp are fixed by assuming the additive quark model relation 
fpN — J-kN together with the Donnachie-Landshoff fit to ttN scattering data. The 
remaining parameters <j p n and e are then determined by fitting to the p-photoproduction 
datafj. First, however, we summarise the earlier SVD results and the data available. 

We are only interested in data at high enough energies for the phase parameter 77 to be 
small and reasonably well estimated by conventional Regge pole ideas. The most precise 
data on complex nuclei at such energies was obtained by a DESY-MIT group ]15|| at a 



photon energy of 6.6 GeV and by a Cornell group pj) at 6.1, 6.5 and 8.8 GeV. The latter 



group also presented results on protons and deuterium, giving measurements of the single 
nucleon cross-section from the same experiment. They then carried out an SVD analysis of 
their own nucleon and nuclear cross-sections, assuming input phases rj which are somewhat 
larger than those given by (|29|), and treating both the 7 — p coupling and the pN total 
cross-section as parameters to be fitted to the data. Their results are summarised in Table 



I. This analysis was repeated by Kroker [24] with phases close to those assumed here, and 
combining both the Cornell and DESY-MIT data at 6.5/6.6 GeV. The resulting parameter 
values are given in Table II, and the corresponding prediction is compared with the data at 
6.1 GeV in Figure ||]. Since both analyses used exactly the same data at 6.1 and 8.8 GeV, 
a comparison of Tables I and II at these energies explicitly confirms the insensitivity of the 
results to small changes in the input phases. 

Here, we repeat this anaysis using the GVD model described above, using the experi- 
mental 7 — p coupling value (^). Since this involves non-trivial numerical computation, we 



4 The contrasting treatment of the cross-section a p N and phase r\ is justified because the results 
are relatively insensitive to small changes in the phase, which is about rj = 0.2 at the energies of 
the data used, but very sensitive to ct p n- 
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simplify the determination of the two free parameters e and <j p n by requiring that the GVD 
predictions exactly reproduce the succesful SVD results on single nucleons; and then deter- 
mine the remaining parameter by a fit to the nuclear data. Satisfactory fits are obtained 
in this way at all three energies 6.1,6.5 and 8.8 GeV as shown in Figures |6], [7| and |8|. The 
corresponding parameter values are shown in Table III. As can be seen, the values obtained 
for the off-diagonal parameter e at the three energies are consistent with each other and with 
the value (||) required by the successful treatment of nuclear shadowing in real photoabsorp- 
tion and deep inelastic scattering ||. The central values for <j p n are slightly larger than 
those obtained in the SVD model fits(cf. Table II), but are consistent within the quoted 
uncertainties arising from the errors on the experimental data. They are also slightly larger 
than, but in qualitative agreement with, the crude prediction (|3]) of the aditive quark model, 
which gives a p ^ ~ 27, 27, 26 mb at E n = 6.1, 6.5, 8.8 GeV respectively. 

In short, the model is in good agreement with the data on p-photoproduction on both 
nucleons and nuclei, with a 7 — p coupling consistent with that measured in electron-positron 
annihihation. The form and magnitude of the correction to simple vector dominance, char- 
acterized by the parameter e, is consistent with that required by the succesful description 
of shadowing in both real photoabsorption and deep inelastic scattering. 
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TABLE CAPTIONS 



Table I Parameters of the SVD fits of the Cornell group [16] to their own data, together 
with the input values of the phase ratio rj. 

Table II Parameters of the SVD fits of Kroker [24] to the Cornell data [16] at 6.1 and 
8.8 GeV, and the combined Cornell and DESY-MIT data [15] at 6.5 GeV, together with the 
input values of the phase ratio rj. 

Table III Parameters of our optical model GVD fits to the Cornell data [16] at 6.1 and 
8.8 GeV, and the combined Cornell and DESY-MIT data [15] at 6.5 GeV. The phase ratio 
i] is given by (29) and the 7 — p coupling is fixed at its measured value (5). 
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FIGURE CAPTIONS 



Figure 1 The SVD diagram for p° photoproduction(/e/t), together with the GVD cor- 
rections to it arising from diffraction dissociation terms(right) . 

Figure 2 Finite width corrections to SVD associated with the 2n- channel. The first 
diagram is a propagator correction, the second a 2it- scattering contribution. 

Figure 3 The reduced wave function $ p calculated at zero impact parameter for lead 
nuclei at 6.1 GeV. The parameters are those of our final fit(see Table III below.) 

Figure 4 The profile function T p calculated for lead nuclei at 6.1 GeV. The parameters 
are those of our final fit (see Table III below). 

Figure 5 Comparison between different theoretical predictions and the complex nuclei 
data for E^=Q.l GeV: the SVD-prediction assuming the parameter values of Table II (dashed 
dotted line); the GVD solution in the approximation (28) using the parameter values given 
by Kroker [24] (dashed line); and the full GVD prediction using the same parameters (solid 
line). 

Figure 6 The optical moddel GVD fit to the complex nuclei data at 6.1 GeV, corre- 
sponding to the parameters of Table III. (The deuteron data point shown is not included in 
the fit for obvious reasons.) 

Figure 7 The optical model GVD fit to the complex nuclei data at 6.5 GeV, correspond- 
ing to the parameters of Table III. (The deuteron data point shown is not included in the fit 
for obvious reasons.) 

Figure 8 The optical model GVD fit to the complex nuclei data at 8.8 GeV, correspond- 
ing to the parameters of Table III. (The deuteron data point shown is not included in the fit 
for obvious reasons.) 
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TABLES 



£ 7 [MeV] 


a v 


a pN [mb] 




6.1 


-0.27 




27.5L1.1 




2.48±0.16 


6.5 


-0.27 




27.9±1.3 




2.60±0.20 


8.8 


-0.24 




25.9±1.0 




2.52±0.16 


TABLE I. Parameters of the SVD fits of the Cornell group [16] to their own 


data, together 


with the input values of the phase ratio 77. 








E y [MeV] 


V 


a pN [mb] 




6.1 


-0.24 




27.5±1.1 




2.44±0.16 


6.5 


-0.23 




26.7±2.0 




2.28±0.10 


8.8 


-0.20 




26.2±1.0 




2.52±0.16 



TABLE II. Parameters of the SVD fits of Kroker [24] to the Cornell data [16] at 6.1 and 8.8 
GeV, and the combined Cornell and DESY-MIT data [15] at 6.5 GeV, together with the input 
values of the phase ratio 77. 



E 1 [MeV] 


V 


a pN [mb] 




€ 


6.1 


-0.26 


28.1L1.1 


2.01 


0.19±0.050 


6.5 


-0.24 


28.4±2.0 


2.01 


0.21±0.035 


8.8 


-0.20 


27.9±1.0 


2.01 


0.28±0.046 



TABLE III. Parameters of our optical model GVD fits to the Cornell data [16] at 6.1 and 8.8 
GeV, and the combined Cornell and DESY-MIT data [15] at 6.5 GeV. The phase ratio 77 is given 
by (29) and the 7 — p coupling is fixed at its measured value (5). 



17 



FIGURES 




FIG. 1. The SVD diagram for p° photoproduction( Ze/t), together with the GVD corrections to 
it arising from diffraction dissociation terms (right). 




FIG. 2. Finite width corrections to SVD associated with the 2n- channel. The first diagram is 
a propagator correction, the second a 2tt- scattering contribution. 
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FIG. 3. The reduced wave function <£ p calculated at zero impact parameter for lead nuclei at 
6.1 GeV. The parameters are those of our final fit (see Table III below.) 
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FIG. 4. The profile function T p calculated for lead nuclei at 6.1 GeV. The parameters are those 
of our final fit (see Table III below). 
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FIG. 5. Comparison between different theoretical predictions and the complex nuclei data for 
.E 7 =6.1 GeV: the SVD-prediction assuming the parameter values of Table II (dashed dotted line); 
the GVD solution in the approximation (28) using the parameter values given by Kroker [24] 
(dashed line); and the full GVD prediction using the same parameters (solid line) . 
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FIG. 6. The optical moddel GVD fit to the complex nuclei data at 6.1 GeV, corresponding to 
the parameters of Table III. (The deuteron data point shown is not included in the fit for obvious 



reasons. 
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FIG. 7. The optical model GVD fit to the complex nuclei data at 6.5 GeV, corresponding to 
the parameters of Table III. (The deuteron data point shown is not included in the fit for obvious 
reasons.) 
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FIG. 8. The optical model GVD fit to the complex nuclei data at 8.8 GeV, corresponding to 
the parameters of Table III. (The deuteron data point shown is not included in the fit for obvious 



reasons. 
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